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Nonlinearity and sharp transmission spectra of random ID nonlinear layered structures are com- 
bined together to produce photon pairs with extremely narrow spectral bandwidths. Indistinguish- 
able photons in a pair are nearly unentangled. Also two-photon states with coincident frequencies 
can be conveniently generated in these structures if photon pairs generated into a certain range of 
emission angles are superposed. If two photons are emitted into two different resonant peaks, the 
ratio of their spectral bandwidths may differ considerably from one and two photons remain nearly 
unentangled. 

PACS numbers: 42.50.Dv,42.50.Ex,42.25.Dd 



I. INTRODUCTION 



After the first successful experiment generating photon 
pairs in the nonlinear process of spontaneous parametric 
down-conversion in a nonlinear crystal and demonstra- 
tion of its unusual temporal and spectral properties more 
than 30 years ago J,], properties of the emitted photon 
pairs have been addressed in detail in numerous investi- 
gations 0. The key tool in the theory of photon pairs 
has become a two-photon spectral amplitude [1, 0, S 
that describes a photon pair in its full complexity. At 
the beginning the effort has been concentrated on two- 
photon entangled states with anti-correlated signal- and 
idler-field frequencies that occur in usual nonlinear crys- 
tals. New sources able to generate high photon fluxes 
have been discovered using, e.g., periodically-poled non- 
linear crystals [71, [8|| , four- wave mixing in nonlinear struc- 
tured fibers llOj]. nonlinear planar waveguides [ll] or 
nonlinearities in cavities [l^ . Chirped periodically-poled 
crystals opened the door for the generation of signal 
and idler fields with extremely wide spectral bandwidths 
[T^ . [l^ . Photon pairs with wide spectral bandwidths 
can also be generated in specific non-coUinear geometries 
p^ . Very sharp temporal features are typical for such 
states that can be successfully applied in metrology (see 
p^6j for quantum optical coherence tomography). Later, 
even two-photon states with coincident frequencies have 
been revealed. They can be emitted provided that the 
extended phase matching conditions (phase matching of 
the wave-vectors and group-velocity matching) are ful- 
filled [l^, [3| • Also other suitable geometries for the 
generation of states with coincident frequencies have been 
found. Nonlinear phase matching for different frequen- 
cies at different angles of propagation of the interacting 
fields can be conveniently used in this case [13, Hi], . 
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Or wave-guiding structures with transverse pumping and 
counter-propa g atin g sig nal and idler fields can be ex- 
ploited ll, [11 ia ii, [11, ii, i^. The last two 

approaches are quite fiexible and allow the generation of 
states having two-photon amplitudes with an arbitrary 
orientation of main axes of their (approximately) gaus- 
sian shape. These approaches emphasize the role of a 
two-photon state at the boundary between the two men- 
tioned cases. This state is, a bit surprisingly, unentangled 
in frequencies and, at the first si ght, should not be very 
interesting. The opposed is true |3l[ . [3^ due to the fact 
that the signal and idler photons are completely indistin- 
guishable and perfectly synchronized in time for pulsed 
pumping. Non-coUinear configurations of bulk crystals of 
a suitable length and pump beam with a suitable waist 
can also be used to generate such state [H, [sl] . These 
states with completely indistinguishable photons are very 
useful in many quantum-information protocols (e.g., in 
linear quantum computing) that rely on polarization en- 
tanglement. 

It has been shown that spectral entanglement af- 
fects polarization entanglement and causes the lowering 
of polarization-entanglement performance in quantum- 
information protocols. For example, the role of spectral 
entanglement in polarization quantum teleportation has 
been analyzed in detail in [sj. We note that this anal- 
ysis is also appropriate for quantum repeaters, quantum 
relays or quantum repeaters. If two photons are identical 
and without any mutual entanglement, the polarization 
degrees of freedom are completely separated from the 
spectral ones and in principle the best possible perfor- 
mance of quantum information protocols is guaranteed. 
In practice, mode mismatches both in spectral and spa- 
tial domains may occur and degrade the performance. It 
has been shown in (35| that Gaussian distributed photons 
with large bandwidths represent states with the best tol- 
erance against errors. Under these conditions high visi- 
bilities of interference patterns created by a simultaneous 
detection of in principal many photons are expected. 

As shown in this paper, nonlinear layered struc- 
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tures with randomly positioned boundaries are a natu- 
ral source of unentangled photon pairs. If the down- 
converted photons are generated under identical condi- 
tions (i.e. into the same transmission peaks) they are in- 
distinguishable and thus ideal for quantum-information 
processing with polarization degrees of freedom. More- 
over their spectra are very narrow and temporal wave 
packets quite long and so their use in experimental se- 
tups is fault tolerant against unwanted temporal shifts 
|3a |. Also because the down-converted photons are gen- 
erated into localized states with high values of electric- 
field amplitudes higher photon-pair generation rates are 
expected. This is important for protocols exploiting sev- 
eral photon pairs. These states with very narrow spectra 
can also be useful in entangled two-photon spectroscopy. 
Superposition of photon pairs generated under different 
emission angles is possible and two-photon states coinci- 
dent in frequencies can be engineered this way. 

We note that the use of high electric-field amplitudes 
occurring in localized states in random structures for the 
enhancement of nonlinear processes has been studied in 
detail for second- harmonic generation in [H, [s^l . 

The paper is organized as follows. The model of spon- 
taneous parametric down-conversion and formulas pro- 
viding main physical characteristics of photon pairs are 
presented in Sec. II. Properties of random ID layered 
structures important for photon-pair emission are stud- 
ied in Sec. III. A typical random structure as a source of 
photon pairs is analyzed in Sec. IV. A scheme providing 
two-photon states coincident in frequencies is analyzed 
in Sec. V. Emission of a photon pair non-degenerate in 
frequencies is investigated in Sec. VI. Sec. VII contains 
Conclusions. 



II. SPONTANEOUS PARAMETRIC 
DOWN-CONVERSION AND 
CHARACTERISTICS OF PHOTON PAIRS 



Spontaneous parametric down-conversion in layered 
media can be described by a nonlinear Hamiltonian H 
in the formalism of quantum mechanics; 



H{t) = eo dr 
Jv 

d(r) : [E(+)(r,t)E(-)(r,i)E(-)(r,i) +h.c.] . (1) 

Vector properties of the nonlinear interaction are char- 
acterized by a third-order tensor of nonlinear coefficients 
d [111. Symbol Ep^^ denotes a positive-frequency pump- 
field electric-field amplitude, whereas negative-frequency 
electric-field amplitude operators eL for m = s,i de- 
scribe the signal and idler fields. Symbol eq is permit- 
tivity of vacuum, V interaction volume, h.c. stands for 
a hermitian conjugated term, and : means tensor multi- 
plication with respect to three indices. The electric-field 
amplitudes of the pump, signal, and idler fields can be 



conveniently decomposed into forward- and backward- 
propagating monochromatic plane waves when describ- 
ing scattering of the considered fields at boundaries of a 
layered structure [H, HO, |4l|. Vector properties of the 
electric-field amplitudes are taken into account using de- 
composition into TE and TM modes. A detailed theory 
has been developed in [s^ and this theory is used in our 
calculations. We note that a suitable choice of vector 
and photonic properties can result in the generation of 
photon pairs with a two-photon spectral amplitude anti- 
symmetric in the signal- and idler- field frequencies [4^ . 
Some of the properties of photon pairs generated in lay- 
ered structures are common with these of photon pairs 
originating in nonlinear crystal super-lattices composed 
of several nonlinear homogeneous pieces [i^, |4j . 

In further considerations, we assume fixed polariza- 
tions of the signal and idler fields and restrict our- 
selves to a scalar model. Then the perturbation solution 
of Schrodinger equation for the signal- and idler-fields 
wave-function results in the following two-photon wave- 
function 0,13]: 

/•oo /-oo 



X dl{u!s)al{uji) exp{iujst) exp(iujit)\vac) . (2) 

The two-photon spectral amplitude (j){u!s,uji) giving a 
probability amplitude of generating a signal photon at 
frequency LUg and its idler twin at frequency LUi is deter- 
mined for given angles of emission and polarizations of 
the signal and idler photons. Creation operator al^{ujrn) 
adds one photon at frequency ujm into field m [i = s,i) 
and |vac) means the vacuum state for the signal and idler 
fields. Both photons can propagate either forward or 
backward outside the structure as a consequence of scat- 
tering inside the structure. Here we analyze in detail the 
case when both photons propagate forward and note that 
similar behavior is found in the remaining three cases. 
More details can be found in [H, HI] . 

The two-photon spectral amplitude <f> can be decom- 
posed into the Schmidt dual basis with base functions 
fs.n and fi^n [30, 45, 46] in order to reveal correlations 
(entanglement) between the signal and idler fields: 



(3) 



A„ being coefficients of the decomposition. Entropy of 
entanglement S defined as [isj 



S 



log2 \ 



(4) 



is a suitable quantitative measure of spectral entangle- 
ment between the signal and idler fields. Symbol log2 
stands for logarithm with base 2. An additional measure 
of spectral entanglement can be expressed in terms of the 
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cooperativity parameter K [46 



K 



1 



(5) 



The two-photon spectral ampUtude (f>{uJs, ^i) describes 
completely properties of photon pairs. Number N of the 
generated photon pairs as well as signal-field energy spec- 
trum Ss{li->s) can simply be determined using the two- 
photon spectral amplitude (soj : 



N 



duJs / duJi\(j){ujs,uJi)\ 
Jo 



duji\(t>{u0s,uii)\ 



(6) 
(7) 



Fourier transform (l){ts,ti) of the spectral two-photon 
amplitude (l){ujs,'-^i), 



2^ I'OD f'OO 

Htsyti) ^ ^ J dojs J duj. 



UJsOJj 



X exp(— zcjs^s) exp(— iwi^i), 



(8) 



is useful in determining temporal properties of photon 
pairs. Symbol u!^ in Eq. ([5]) denotes the central fre- 
quency of field TO, m = s,i. This Fourier transform is 
linearly proportional to the two-photon temporal ampli- 
tude A defined along the expression 

A{ts,U) = {vac\Ei+'>{0,to + ts)^+\0,to + U)\^/Js,^{tn)) 

(9) 

and giving a probability amplitude of detecting a signal 
photon at time together with its twin at time U. We 
have [Sg: 



AnenCD 



(10) 



where c is speed of light in vacuum and B transverse 
area of the fields. Photon fiux Af of, e.g., the signal field 
is determined using a simple formula (valid for narrow 
spectra) (soj : 



(11) 



Direct measurement of temporal properties of photon 
pairs is impossible because of short time scales needed 
and so temporal properties may be detected only indi- 
rectly. Time duration of the (pulsed) down-converted 
fields as well as entanglement time can be experimen- 
tally addressed using Hong-Ou-Mandel interferometer. 
Normalized coincidence-count rate B^'~'^ in this inter- 
ferometer is derived as follows H: 



i?r^(ro 

pin) 



1 - p{n), 
1 



(12) 



i?(0,0) 



Re 



duj^ 



dlOi UJsLUi 



{ujs,uj.i)4}*{uji,ujs) exp^iuji Ti ) exp {-iujsTi) 



(13) 



Normalization constant i?(0, 0) occurring in Eq. (fT3|) is 
given in Eq. (fT5|l bellow. Relative time delay ti between 
the signal and idler photons changes in the interferome- 
ter. 

A preferred direction of correlations between the 
signal- and idler-field frequencies can be determined from 
the orientation of coincidence-count interference fringes 
in Franson interferometer [25| that is characterized by 
normalized coincidence-count rate in the form: 



1 



1 



Re {2R{ts,0) + 2R{0,n) 

(14) 



4 8i?(0,0) 

+R{Ts,n) + R{Ts,-n)}. 

The function R in Eq. is defined as follows 

R{Ts,n) ^ 



dhJs I duji uJsUJi\(l){uJs,uJi)\ 
Jo 

X exp(iWsTs) exp(iLi;iTi). (15) 

Time delay (m — s, i) corresponds to a relative phase 
shift between two arms in the path of photon to. 



III. PROPERTIES OF RANDOM ID LAYERED 
STRUCTURES 

We consider a ID layered structure composed of two 
dielectrics with mean layer optical thicknesses equal to 
Ao/4 (Ao = 1 X 10~^ m is chosen). Such structure, for 
example, can be fabricated by etching a crystal made 
of LiNbOs and filling free spaces with a suitable ma- 
terial (Si02). The optical axis of LiNbOs is parallel 
to the planes of boundaries and coincides with the di- 
rection of fields' polarizations that correspond to TE 
modes. This geometry exploits the largest nonlinear 
coefficient of tensor d (dzzz)- The fields propagate as 

ordinary waves with the dispersion formula ri^{uj) = 

4.9048-1-0.11768/(354.8084-0.0475^2) - 9.6398^2 
The second material (Si02) is characterized by its index 
of refraction n(Lo) = 1.45. A random structure is gener- 
ated along the following algorithm: 

1. The number A'cicm of elementary layers of optical 
thicknesses Ao/4 is fixed. Then the material of each 
elementary layer is randomly chosen. 

2. At each boundary between two materials, an ad- 
ditional random shift of the boundary position 
governed by a gaussian distribution with variance 
Ao/40 is introduced. 

The randomness given by a random choice of the material 
of each layer is crucial for the observation of features typ- 
ical for random structures. The randomness introduced 
by gaussian shifts of boundary positions is only additional 
and does not modify considerably properties of a random 
structure. It also describes errors occurring in the fabri- 
cation process. We note that the same generation algo- 
rithm has been used in (36l . ISTj where second-harmonic 
generation in random structures has been addressed. 
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The process of spontaneous parametric down- 
conversion can be efficient in these structures pro- 
vided that the signal and idler fields are generated into 
transmission peaks (corresponding to spatially localized 
states) and the structure is at least partially transparent 
for the pump field. This imposes restrictions to the al- 
lowed optical lengths of suitable structures. They have to 
have such lengths that the down-converted fields (usually 
degenerate in frequencies) are localized whereas there oc- 
curs no localization at the pump-field wavelength. This 
is possible because the shorter the wavelength the larger 
the localization length [i^ for a given structure. We as- 
sume the pump-field (signal-, idler-field) wavelength Xp 
(As, Xi) in the vicinity of Ao/2 (Aq). 

Numerical simulations for the considered materials and 
wavelengths have revealed that the best suitable numbers 
-^cicm of elementary layers lie between 200 and 400. In 
general, the greater the contrast of indices of refraction 
of two materials the smaller the number of needed layers. 
Widths of transmission peaks corresponding to localized 
states in a random structure may vary by several orders of 
magnitude. As an example, probability distribution Pax^ 
of widths of transmission peaks for an ensemble of struc- 
tures with A'cicm = 250 elementary layers (optical lengths 
of these structures lie around 6 x 10 ~^ m) is shown in 
Fig. [T^. Comparison with the distributions appropriate 
for iVoicm — 500 and Ndcm = 750 in Fig. [1^ demonstrates 
that the longer the structure the narrower transmission 
peaks can be expected. Localization optical length at 
Ao [13, nil in the direction perpendicular to the bound- 
aries lies around 22 x 10^*^ m for the considered structure 
lengths. As for the simulation, we have randomly gener- 
ated several tens of thousands structures in order to have 
around 3 x lO** transmission peaks for the statistical anal- 
ysis. Transmission peaks can be easily distinguished from 
their flat surroundings for which intensities lower than 
one percent of the maximum peak intensity are typical. 
Widths of intensity transmission peaks have been deter- 
mined as FWHM. Simulations have revealed that peaks 
with very small intensity transmissions prevail, but there 
occur also peaks with intensity transmissions close to one. 
Such peaks are useful for the generation of photon pairs 
because of the presence of high electric-field amplitudes 
inside the structure. Photon pairs can be also generated 
under nonzero angles of emission. The greater the radial 
angle 9s of signal-photon emission the narrower the trans- 
mission peaks as documented in Fig. [TJ). Also the local- 
ization lengths (projected into the direction perpendicu- 
lar to the boundaries) shorten with an increasing radial 
angle Og of signal-photon emission. Localization length 
equal to 9.5 x 10~^ m (2.5 x 10~^ m) has been deter- 
mined at the angle of signal-photon emission 9s = 30 deg 
(60 deg). 

Fabrication of such structures is relatively easy due 
to allowed high tolerances. A typical sample has several 
transmission peaks at different angles 9^ of emission for a 
given frequency of the signal field. If we assume that also 
the idler field is tuned into the same transmission peak 



a) 15 -| 
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AAs (10 m) 

FIG. 1: Probability distribution PaXs of having a transmis- 
sion peak with width (FWHM) AAs in the signal field at cen- 
tral wavelength near Ao (a) for an ensemble with A'cicm = 250 
(solid line), A'eiem = 500 (solid line with *), and A^dcm ~ 750 
(solid line with A) at 6s = deg and (b) for angles 6s of 
signal-photon emission deg (solid line), 30 deg (solid line 
with *), and 60 deg (solid line with A) using an ensemble 
with N^icn-i = 250. 



the (normally incident) pump field has to have twice the 
frequency corresponding to this peak. In non-coUinear 
geometry the transverse wave- vectors of the signal and 
idler fields have to have the same magnitudes and op- 
posed signs in order to fulfill phase-matching conditions 
in the transverse plane. If the angle 9s of signal-photon 
emission increases a given transmission peak survives in- 
creasing its central frequency cj^. Moreover the depen- 
dence of central frequency a;^ on radial emission angle 
9s can be considered to be linear in a certain interval 
of angles 9s- This property leads to wide tunability in 
frequencies. 



IV. PROPERTIES OF THE EMITTED PHOTON 
PAIRS 

We assume a normally incident pump beam forming a 
TE wave and generation of the TE-polarized signal and 
idler fields into the same transmission peak. Thus both 
emitted photons have the same central frequencies and 
the central frequency of the pump field is twice this fre- 
quency, i.e. the conservation law of energy is fulfilled. If a 
structure is pumped by a broadband femtosecond pulse, 
photon pairs are generated into a certain range of emis- 
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FIG. 2; Contour plot of relative energy spectrum S^^ of the 
signal field versus normalized signal-field frequency 2aja/a;p 
and angle Qs of signal-field emission. The relative spectrum 
S^f^ is given by the ratio of actual spectrum of the structure 
and spectrum of a reference homogeneous structure contain- 
ing the same amount of perfectly-phase-matched nonlinear 
material. The structure is pumped by a pulse 250 fs long. 



sion angles. If a signal photon occurs at a given radial an- 
gle Qs its idler twin is emitted at its radial angle Oi — —Og 
and wave- vectors of both photons and the pump field lie 
in the same plane (i.e., their azimuthal angles '0s and ipi 
coincide) as a consequence of phase matching conditions 
in the transverse plane. The central frequency uj^ of a 
signal photon depends on the angle Og of signal-photon 
emission, as illustrated in Fig. [5] showing the signal- field 
intensity spectrum 5™^. Width of the spectrum Sg coin- 
cides with the width of intensity transmission peak be- 
cause the transmission peaks are very narrow and so all 
frequencies inside them have nearly the same conditions 
for the nonlinear process. This means that linear proper- 
ties of the photonic-band-gap structure have a dominant 
role in the determination of spatial properties of the gen- 
erated photon pairs. Localization of the signal and idler 
fields inside the considered photonic-band-gap structure 
leads to the enhancement of photon-pair generation rate 
up to 5000 times (in the middle of transmission peak) as 
can be deduced from Fig. [21 This is in accordance with 
the finding that the process of second-harmonic genera- 
tion can be enhanced by 3 or 4 orders of magnitude [s^l 
in these structures. This enhancement is at the expense 
of dramatic narrowing of the range of the allowed fre- 
quencies of emitted photons. We estimate from 10 to 
lO'^ generated photon pairs into the whole emission cone 
per 100 mW of the pump power depending on the width 
of transmission peak. The wider the transmission peak 
the higher number of pairs is expected. We note that 
extremely narrow down-converted fields can be obtained 
also from parametric pumping of ^''Rb atoms [soj . 

Two-photon spectral amplitudes ^(ws, Wj) for different 
angles Og of signal-photon emission are very similar in 
their shape; they differ in their central frequencies. We 
note that the central frequencies of the signal and idler 



FIG. 3; Contour plot of squared modulus \4>'^ of the two- 
photon spectral amplitude for the transmission peak at 6s — 
deg. The structure is pumped by a pulse 250 fs long; 101^ is 
normalized such that 4 J dujs J dLJi\4>{ujs,uJi)\^ /(cdp)^ = 1. 
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FIG. 4: Squared modulus \fn\^ of mode functions with n=l 
(solid line), 2 (solid line with *), and 3 (solid line with A) of 
Schmidt decomposition of the two-photon spectral amplitude 
in Fig.O /„,s(lj) — f„,i{cL>) = fn{i-o) due to symmetry. Mode 
functions /„ are normalized such that 2 J dtt;|/„(tij) p/o;" — 1. 



photons are the same {uj'^ — uj^) for a given angle 9g 
of emission. A typical shape of the two-photon spec- 
tral amplitude 4'{ujg, LUi) resembling a cross in its contour 
plot (see Fig. [3]) reflects the fact that the signal and idler 
fields are nearly perfectly separable. This is confirmed 
by Schmidt decomposition of the amplitude 4>{ujg,uji) in 
which only the first mode is important (Ai = 1.00). En- 
tropy S of entanglement for the two-photon spectral am- 
plitude 4> in Fig. [3] equals 0.00 and cooperativity param- 
eter K is 1.00. The spectral dependence of the first three 
mode functions in the decomposition is shown in Fig. IH 
A typical temporal two-photon amplitude A{tg,ti) 
spreads over tens or hundreds of ps reflecting narrow fre- 
quency spectra of the signal and idler fields. Its contour 
plot resembles a droplet ^] that originates in the zig-zag 
movement of the emitted photons inside the structure 
that delays the occurrence time of photons at the out- 
put plane of the structure. Both photons have the same 
wave-packets due to identical emission conditions as can 
be verified in Hong-Ou-Mandel interferometer showing 
the visibility equal to one. 
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FIG. 5: Spread of the correlation area in radial {aA6i, solid 
line) and azimuthal {aA^pi, solid line with *) directions of 
idler-photon emission as a function of pump-beam diameter a. 
The transmission peak in the structure with 250 elementary 
layers occurs at S° = 26.3 deg, i.e. the central radial idler- 
photon emission angle equals - 26.3 deg. Logarithmic scale 
on the X axis is used. 



Assuming the normally incident pump beam the 
signal-field intensity profile at a given frequency in the 
transverse plane is nonzero around a circle due to the 
rotational symmetry of the photonic-band-gap structure 
around the pump-beam propagation direction. This sym- 
metry assumes an appropriate rotation of polarization 
base vectors. However, values of intensity change around 
the circle depending on fields' polarizations (defined by 
analyzers in front of detectors) and properties of nonlin- 
ear tensor d. 

Correlations between the signal and idler fields in the 
transverse plane are characterized by a correlation area 
that is defined by the probability of emitting an idler 
photon in a given direction [defined by angular declina- 
tions AOi (radial direction) and Atpi (azimuthal direc- 
tion) from the ideal direction of emission] provided that 
the direction of signal-photon emission is fixed. Correla- 
tion area is in general an ellipse with typical lengths in 
radial (ctA^^) and azimuthal (aA-tpi) directions. These 
angles are very small for plane- wave pumping, typically 
of the order of 10^'^ — 10^^ rad. However, focusing of 
the pump beam can increase their values considerably 
as shown in Fig. [51 where the pump-beam diameter a 
varies from 30 fj,m up to 1 mm. Spread of the correla- 
tion area in radial direction is smaller compared to that 
in azimuthal direction, because emission of a photon is 
more restricted in radial direction (for geometry reasons) 
by narrow bands of the photonic structure. Release of 
strict phase-matching conditions in the transverse plane 
caused by a focused pump beam (with a circular spot) af- 
fects radial and azimuthal angles of emission in the same 
way. 

We note that spreading of the signal- and idler-field in- 
tensity profiles caused by pump-beam focusing has been 
experimentally observed in fsij for a type-II bulk crystal. 



V. TWO-PHOTON STATES COINCIDENT IN 
FREQUENCIES 

Superposition of photon pairs with signal photons 
emitted under different radial emission angles 9s results 
in states coincident in frequencies. This spectral beam 
combining of fields from different spatial modes and with 
different spectral compositions can be achieved using an 
optical dispersion element like a grating [s^]- This tech- 
nique has already been applied in the construction of a 
source of photon pairs that uses achromatic phase match- 
ing and spatial decomposition of the pump beam |2l[ . 

There is approximately a linear dependence between 
the central frequencies tj^ — uj^ and radial angle 9s of 
signal-field emission for sufficiently large values of dg as- 
suming a normally incident pump beam. The resultant 
two-photon spectral amplitude ^M{,0Js,uJi) after super- 
posing photon pairs from M equidistantly positioned pin- 
holes present both in the signal and idler beams can be 
approximately expressed as: 

M-l 

$A/('^s,'^i) — e'xjp{i(pn)(l){LUs + nAuj,LUi + nAuj), 

Tl = 

(16) 

where Alu is the difference between the signal-field cen- 
tral frequencies of fields originating in adjacent pinholes. 
Phase ip determines the difference in phases of two two- 
photon spectral amplitudes at the central signal- and 
idler-field frequencies coming from adjacent pinholes. 
The two-photon spectral amplitude (l){ujs,uji) in Eq. (|16p 
belongs to a two-photon state coming from the first pin- 
hole. Schmidt decomposition of the two-photon spec- 
tral amplitude describing photon pairs coming from 
M pinholes shows that there are nearly M independent 
modes [determined by the value of cooperativity param- 
eter K introduced in Eq. ([5])]. These modes are collec- 
tive modes, i.e. they have non-negligible values in areas 
of frequencies associated with every pinhole (see Figs. [H] 
and [7] for M — 2 pinholes) . Such states are perspective 
for the implementation of various quantum-information 
protocols. 

The Fourier transform $j\/(is, ti) defined in Eq. ([5]) can 
be approximately expressed as follows: 

$j\/(i.,t»)« ^ exp[i{<fi-Auj{ts+t,))n](f>{ts,U), (17) 

rt=0 

where 

(t){ts,U) = — dujs / dwi(/)(ws,Wi)exp(iWsts) 

X e-x:p[iujiti) (18) 

means the Fourier transform of a two-photon ampli- 
tude associated with one pinhole. The expression in 
Eq. (|17p indicates that interference of photon-pair ampli- 
tudes originating in different pinholes creates interference 
fringes in the sum ts -\- ti of the occurrence times of signal 
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FIG. 6: Contour plot of squared modulus |"I>i\/|^ of the two- 
photon spectral amplitude created by superposing photon- 
pair amplitudes from M — 2 pinholes. The structure is 
pumped by a pulse 250 fs long; |"1?mP is normalized such 
that AjduJsJ diJi\<^M{ujs,LOi)\^/(tOpf = 1. 




FIG. 7: Squared modulus |/np of mode functions with n=l 
(solid line), 2 (solid line with *), and 3 (solid line with A) 
of Schmidt decomposition of the two-photon spectral am- 
plitude $A/ in Fig. [6l fn,a{i^) = fn,i{i^) = fn{tLl) duC tO 

symmetry. Mode functions /„ are normalized such that 

2 / dLv\U{L0)\^/LV° = 1. 



and idler photons. This is caused by positive correlations 
in the signal- and idler-field frequencies. The analysis of 
the sum in Eq. (|17p shows that the greater the number 
M of pinholes, the narrower the range of allowed val- 
ues of the sum tg + U [lim7\/^oo 'Ylin=o exp(«x7i) « 6{x)\. 
Comparison of shapes of the two-photon temporal am- 
plitudes AM{ts,ti) derived in Eq. pO)) for M — 2 and 
Af = 8 pinholes shown in Fig. [5] reveals this tendency 
for localization in time domain. These features originat- 
ing in positive correlations of the signal- and idler-field 
frequencies are experimentally accessible by measuring 
the pattern of coincidence-count rate i?^ in Franson in- 
terferometer for sufficiently large values of signal- and 
idler-photon delays and r^. A (nearly separable) two- 
photon state from one pinhole creates a chessboard tilted 
by 45 degrees (see Fig. Oi). If amplitudes from several 
pinholes are included typical fringes oriented at 45 de- 
grees become visible in coincidence-count patterns given 
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FIG. 8: Contour plot of squared modulus P of the tempo- 
ral two-photon amplitude for (a) M = 2 and (b) M = 8. The 
corresponding two-photon spectral amplitude $m for M = 2 
is plotted in Fig. [S] Positions of the edge pinholes for Af = 8 
coincide with those for M — 2 and the remaining six pinholes 
are equidistantly distributed in-between them. Normalization 
is such that J dta J dti\AM{ts,ti)\'^ = 1. 



by rate i?^ (compare Figs. and[9j3). The greater the 
number M of pinholes, the better the fringes are formed. 

Superposition of photon-pair amplitudes (with a suit- 
able phase compensation) from a given range of signal- 
field emission angles 9s (and the corresponding range 
of idler-field emission angles 9i) as defined by rectan- 
gular apertures gives a two-photon spectral amplitude 
$ having a cigar shape with coincident frequencies (see 
Fig. nop . The phase compensation requires the intro- 
duction of an additional phase shift depending on the 
radial emission angle 9^ and can be accomplished, e.g., 
by a wedge-shaped prism. Schmidt decomposition of the 
two-photon amplitude $ shows that several modes are 
non-negfigibly populated (Ai — 0.78, A2 = 0.13, and 
A3 = 0.04). Typical profiles of mode functions /„ for this 
configuration are plotted in Fig. [TT] The first and most 
populated mode extends over all frequencies whereas os- 
cillations are characteristic for the other modes. The 
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FIG. 9: Contour plot of normalized coincidence-count rate 
in Franson interferometer as a function of signal- and idler- 
field delays Ts and Ti for (a) M = 1 and (b) M — 8 pinholes. 
Positions of pinholes are defined in caption to Fig. [S] 



larger the spectral width of the signal and idler fields 
the larger the cooperativity parameter K. This means 
that the number of effective independent modes can be 
easily changed just by changing the height of rectangu- 
lar apertures. This makes this source of photon pairs 
extraordinarily useful. Both emitted photons are per- 
fectly indistinguishable providing visibility equal to one 
in Hong-Ou-Mandel interferometer. Positive correlation 
in the signal- and idler-field frequencies gives coincidence- 
count interference fringes in Franson interferometer tilted 
by 45 degrees for sufficiently large values of the signal- 
and idler-field delays and Ti. 



VI. SPECTRALLY NON-DEGENERATE 
EMISSION OF PHOTON PAIRS 

The simplest way for the observation of spectrally 
non-degenerate emission of a photon pair is to consider 
collinear geometry and exploit a random structure with 
two different transmission peaks. Pumping frequency is 



FIG. 10: Contour plot of squared modulus |$|^ of the 
two-photon spectral amplitude created by superposition of 
photon-pair amplitudes from a certain range of signal-field 
(and idler-field) emission angles 6s (Oi) with a suitable 
phase compensation. The structure is pumped by a pulse 
250 fs long; |$|^ is normalized according to the condition 
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FIG. 11: Squared modulus I/71P of mode functions with 
n=l (solid line), 2 (solid line with *), and 3 (solid line 
with A) of Schmidt decomposition of the two-photon spec- 
tral amplitude $ in Fig. [TD] fn,s{'.o) — fn,i{u)) = fn{u>) due 
to symmetry. Mode functions /„ are normalized such that 



2 / da;|/„(w)|Va 



then given as the sum of central frequencies of these 
transmission peaks into which signal and idler photons 
are emitted. The generation of a suitable structure us- 
ing the algorithm and geometry presented in Sec. Ill in 
this case is by an order of magnitude more difficult com- 
pared to that providing photon pairs degenerated in fre- 
quencies. Structures having two peaks with considerably 
different bandwidths of the transmission peaks are espe- 
cially interesting. Spectral bandwidths of the signal and 
idler fields differ accordingly. Such states are interest- 
ing in some applications, e.g., in constructing heralded 
single-photon sources. 

The two-photon spectral amplitude 4>{ujs,uji) has a 
cigar shape prolonged along the frequency of the field 
with a larger spectral bandwidth (say the signal field). 
On the other hand, the two-photon temporal amplitude 
A{ts, ti) is broken into several islands along the axis giv- 
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FIG. 12: Photon fluxes of the signal {Afs, solid line) and idler 
{Mi, solid line with *) fields for a structure with Ndcm = 
250 layers having two different transmission peaks in coUinear 
geometry. The ratio of their intensity bandwidths equals cca 
4. The structure is pumped by a pulse 250 fs long. 
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FIG. 13: Coincidence-count rate in Hong-Ou-Mandel 

interferometer as a function of relative time delay n for the 
two-photon state used in Fig. \T2\ 

ing the idler-photon detection time t^. Phase modula- 
tion of the two-photon spectral ampHtude (p{u>s,uji) with 
faster changes along the frequency LOi and slower changes 
along the frequency is responsible for this behavior. 
Photon fluxes of the signal and idler flelds are then com- 
posed of several peaks as documented in Fig. [T^ This 
feature is reflected in an asymmetric dip in coincidence- 
count rate i?^*-*^ in Hong-Ou-Mandel interferometer, 
that also shows oscillations at the difference of the cen- 



tral signal- and idler-field frequencies. Asymmetry of the 
interference dip shown in Fig. [T^] can be related to the 
shape of spectral two-photon amplitude (j) Origin of 
these effects lies in delays caused by the zig-zag move- 
ment of the generated photons inside the structure. 



VII. CONCLUSIONS 

Nonlinear random layered structures in which an op- 
tical analog of Anderson localization occurs have been 
analyzed as a suitable source of photon pairs with nar- 
row spectral bandwidths and perfect indistinguishability. 
Spectral bandwidths in the range from 1 nm to 0.01 nm 
are available for different realizations of a random struc- 
ture. Photon pairs with the same signal- and idler-field 
bandwidths as well as with considerably different band- 
widths can be generated. Random structures are ffexi- 
ble as for the generated frequencies and emission angles. 
Two-photon states with coincident frequencies and vari- 
able spectral bandwidth can be reached if two-photon 
amplitudes of photon pairs generated into different emis- 
sion angles are superposed. Also two-photon states with 
signal- and idler-field spectra composed of several peaks 
and characterized by collective spectral mode functions 
are available. All these states are very perspective for 
optical implementations of many quantum information 
protocols. Possible implementation of these sources into 
integrated optoelectronic circuits is a great advantage. 
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